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Nonproduct n-qubit pure states with maximum dimensional stabilizer subgroups of the group 
of local unitary transformations are precisely the generalized n-qubit Greenberger-Horne-Zeilinger 
states and their local unitary equivalents, for n > 3, n 7^ 4. We characterize the Lie algebra of the 
stabilizer subgroup for these states. For n — 4, there is an additional maximal stabilizer subalgebra, 
not local unitary equivalent to the former. We give a canonical form for states with this stabilizer 
as well. 



00 

o 
o 

(N 



Oh! 



> 
in 
o 



o 
o 



X 



PACS numbers: 03.67.Mn 



I. INTRODUCTION 



The desire to measure and classify entanglement for 
multiparty states of n-qubit systems has been moti- 
vated by potential applications in quantum computation 
and communication that utilize entanglement as a re- 
source [l|, Q . Entanglement classification is also an in- 
teresting question in the fundamental theory of quantum 
information. 

It is natural to consider two states to have the same 
entanglement type if one is transformable to the other 
by a local unitary transformation [1,0]. Thus any entan- 
glement measure, as a function on state space, must be 
invariant under the action of the local unitary group. One 
such invariant is the isomorphism class of the stabilizer 
subgroup (the set of local unitary transformations that do 
not alter a given state) and its Lie algebra of infinitesimal 
transformations. Inspired by ideas originally laid out for 
3-qubit systems 0, @ , the authors of the present paper 
have achieved a number of results for systems of arbi- 
trary numbers of qubits on the structure of stabilizer Lie 
subalgebras and in particular, those stabilizers that have 
maximum possible dimension. Evidence that maximum 
stabilizer dimension is an interesting property is the fact 
that states that have such stabilizers also turn out to 
maximize other known entanglement measures andplay 
key roles in quantum computational algorithms [g, |7[. 
In we describe a connection between stabilizer struc- 
ture and the question of when a pure state is determined 
by its reduced density matrices. 

States with maximum stabilizer dimension are prod- 
ucts of singlet pairs (with an uncntangled qubit when the 
number of qubits is odd) and their local unitary equiva- 
lents [1, In the present paper, we show an analogous 
result for nonproduct states with maximum stabilizer di- 
mension. These are the generalized ?T.-qubit Greenberger- 
Horne-Zeilinger (GHZ) states 

a|00---0)+/3|ll---l), a,P^O 



and their local unitary equivalents for 3 or more qubits, 
together with an additional class of states for the spe- 
cial case of 4 qubits. In order to prove these results, we 
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extend the stabilizer analysis of our previous work [11 1 
from stabilizers of pure states to stabilizers of arbitrary 
density matrices. 



II. THE MAIN RESULT 

The n-qubit local unitary group is 

G ~ Go X Gi X ■ X Gn 

where Go = U{1) is the group of phases and Gj = SU{2) 
is the 3-dimensional group of single qubit rotations in 
qubit j for 1 < j < n. Its Lie algebra 



is the set of infinitesimal local unitary transformations 
acting on the tangent level, where go is 1-dimensional 
and each gj = su{2), the set of skew-Hermitian matrices 
with trace zero. Given an n-qubit state vector |^/;), let 

Stab^ = {gGG:5lV^) = |^)} 

denote the subgroup of elements in G that stabilize ip, 
and let 

be the corresponding Lie algebra. 

In fllj we show that the maximum possible dimen- 
sion for is n — 1 when tp is not a product state, for 
n > 3. Further, we show that for n 7^ 4, we must have 
dim PjK^, = 1 for all j, where Pj-. g ^ Qj is the natural 
projection. In fact this condition holds for the general- 
ized n-qubit GHZ state 

a|00---0)+/3|ll---l) 

and its LU equivalents. For n = 4, there is an additional 
possible condition for to attain maximum dimension, 
and that is dim PjK^ = 3 for all j and = su(2). The 
main result of the present paper is the converse of these 



2 



previous results. We show that the generahzed n-qubit 
GHZ states (and LU equivalents) are the only nonprod- 
uct states that have maximum stabilizer dimension, ex- 
cept for the 4-qubit case. For those 4-qubit states with 
stabilizer isomorphic to su{2), we give a canonical rep- 
resentative for each LU equivalence class. Here is the 
formal statement. 

Theorem 1. (Classification of states witii maximal 
stabilizer) Let |?/') be a nonproduct n-qubit state with 
stabilizer subalgebra K.^ of maximum possible dimension 
n — 1, for some n > 3. One of two conditions must hold: 

(i) dimPjKii, ~ 1 for 1 < j < n, or 

(ii) n = 4, dimPjii'^, = 3 for I < j < n and K^, = 
su{2). 

If (i) holds, then is LU equivalent to a generalized n- 
qubit GHZ state. If (ii) holds, then is LU equivalent 
to a state of the form 

IV-) = a(|0011) + |1100)) 
+ 6(|1001) + 10110)) 
+ c(|1010)-h 10101)) 

for some complex coefficients a, 6, c satisfying a > 0, 
abc 7^ and a + 6 + c = 0. Furthermore, this is a unique 
representative of the LU equivalence class for of this 
form. 

The proof divides naturally into two sections, one for 
each of the two conditions on stabilizer in the statement 
of the theorem. We consider conditions (i) and (ii) in 
sections Hvl and Iv l respectively. It is convenient to work 
in terms of density matrices and their local unitary sta- 
bilizers; the following preliminary section establishes the 
necessary extension of our previous results to local uni- 
tary action on density matrices. 

III. PRELIMINARY PROPOSITIONS ON 
STABILIZER STRUCTURE 

For the set of ri-qubit pure and mixed density matrices, 
we may omit the phase factor in local unitary operations, 
and take the local unitary group to be 

S'C/(2)" = Gi X • • • X G„ C G 

and its Lie algebra to be 

su(2)" = 01 © • • • e 0„ C 0. 

Given an n-qubit density matrix p, an element g G 
SU{2Y acts on p by 

9- p = gpg'^- 

An element X G .sm(2)" acts on the infinitesimal level by 
X-p^[X,p]^Xp-pX. 



Let 

Sis^hp = {g e SU{2r : gpg^ =p} 

denote the subgroup of elements that stabilize p, and let 

Kp^{X !E su{2Y : [X, p] = 0} 

be the corresponding Lie algebra. We use multi-index 
notation / = (ii, i2, . . . , i„), where each ik is a binary 
digit, to denote labels for the standard computational 
basis for state space. We write I^ to denote the bitwise 
complement of /. 

One expects a natural correspondence between K^p and 
Kp in the case where p = (-01, and indeed there is. 

Proposition 1. Let p — \^) {tp\ be a pure n-qubit density 
matrix. Then Kp = PK^, where -P: ^ 0i © • • • © 0n is 
the natural projection that drops the phase factor. 

Proof. We begin with a simple observation. Let X E g 
and let = X Skew-Hermicity of X gives us 

{ij\X^-{<P\. (1) 

To see that PK^, C Kp, \ct X e K^,, say X = 
{-it,Xi,X2,...,Xn). Let Y = PX, so y IV-) = it\ip). 
Then 

[Y, p]=Y - Y^itp^ itp - 0. 

To get the second equality, apply ([T]) to = it\ip) — 
Y\^). 

Conversely, we show that Kp C PK^. Let Y G Kp, 
and let |(/)) = F 1-0). Using ([T]) again, we have 

= [Y,p]=Y\ij){i^~\ij){^\Y 

= 10) (^1 + IV) (01 = 10) (V'l + (10) (Vl)t 

so 10) is skew-Hermitian. Thus there is some unitary 
U such that C/|0) (i/jjC/^ is diagonal, with pure imagi- 
nary entries. In fact there is only one nonzero diagonal 
entry, say is for some real s, since |0) (0| has rank one. 
Therefore U |0) = isU so |0) = is |V'). Thus Y is the 
projection of {—is, Yi, Y2, . . . , Y^) in K^. This concludes 
the proof. □ 

Let us use the symbol Pj to denote both the projection 
Pj : — > Qj and also Pj-. su(2)" — > 0^. Applying Pj to 
both sides of Kp = PK^ yields the following corollary. 

Corollary 1. Let p = (-01 be a pure n-qubit density 
matrix. Then PjKp ~ PjK^ for 1 < j < n. 

The next proposition establishes that we may work 
with either Kp or K^ to calculate stabilizer dimension. 

Proposition 2. Let p = \ip) {ip\ be a pure n-qubit density 
matrix. Then dim Kp — dim K^ . 
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Proof. Let Xi, . . . , Xr £ be linearly independent. Let 
Yfe = PXj, for 1 < fc < r and suppose that 

for some scalars Ofc. Since Yk ~ Xk — itk for some tk, we 
have 

^^OkXk = i'^aktk- 

Since ^a^Xfe is in K^, so is i^autk, and so both of 
these sums must be zero. Since the X^. are independent, 
all the afe must be zero, and hence the Yk are also linearly 
independent. □ 



We choose the basis 



A = 



i 

-i 



B = 



1 

-1 



c = 



i 

1 



for su{2). We denote by Aj the element in g that has A in 
the jth qubit slot and zeros elsewhere. By slight abuse of 
notation, we use the same symbol for the corresponding 
element in su(2)". Analogous notation applies for Bj 
and Cj, for 1 < j < n. 

We will make repeated use of the following basic cal- 
culation. Let X ^ J2 tkAk be an element of the local 
unitary Lie algebra su(2)". It is straightforward to check 
that 



= = ^C(/,JK,/|/) {J\ 



(2) 



where 



As a consequence, we have the following. 

Proposition 3. Let p be an n-qubit density matrix, let 
X = "^tkAk, and suppose that X £ Kp. Then either 
C(/, J) = or pi^j = for all I, J. 

We record here a proposition proved in [lo| 
(Lemma 3.8) with an alternative proof using Proposi- 
tion [3l 

Proposition 4. (Stabilizer criterion for an unen- 
tangled qubit) Let p ~ lip) (■01 be a pure n-qubit density 
matrix where lip) ~ X^'^-f 1-^)- U £ Kp, then the £th 
qubit is unentangled. 

Proof. Choose any nonzero state coefficient c/'. Apply 
Proposition [3] to X = Ai, L = L' and J for which jg ^ i'^. 
Since (,{!', J) = ±2i ^ 0, wc conclude that cj must be 
zero. □ 

Proposition 5. (Stabilizer criterion for general- 
ized n-qubit GHZ states) Let = c/ |/) be a 

nonproduct n-qubit state vector for some n>3, and let 



V 



Y^tkAk-. ^tfc = 



We have V = if and only if ip is a generalized n-qubit 
GHZ state. 

Proof. In Section V of [ll[ we show that if -0 is a gener- 
alized n-qubit GHZ state, then A',/, = V . Now we prove 
the converse. 

Suppose that I'C.^, = V . Choose a multi-index / such 
that c/ 7^ 0. If J is another multi-index different from / 
in indices K, ~ {ki, . . . ,kp\ and equal to / in some index 
r, let 



A = ^(-i)^^A,-|y^(-i)' 



so that the sum of the coefficients of the Aj is zero, so X is 
in AT^ . Viewing X as an element of Kp and applying ([2]) , 
we have 



C(/, J) = 2i 



P 

{-lytt = 22 ^(-1)2"^'.- ^2ip^{) 

3 = i 



SO pi. j = cjCj = by Proposition [3l so cj = 0. It 
follows that has the form |?/;) = a |/) -f- /3 for 
some nonzero a, (3. Finally, if / does not consist of all 
zeros or all ones, say ik = 0,ie = 1, then let X = A^ — Ai, 
then (^{I, I'^) = 4i 7^ 0, so pi jc would have to be zero, 
but this is impossible. Thus wc conclude that \^j) ~ 
a|00---0)-h^|ll---l). □ 



IV. MAXIMAL STABILIZERS WITH 
1-DIMENSIONAL PROJECTIONS IN EACH 
QUBIT 

In this section we consider the consequences of con- 
dition (i) in the statement of Theorem [T] By virtue of 
Propositions [T] and [2l we may work with I{p in place of 
K^. If dim PjKp ~ 1 for all j, we can apply an LU trans- 
formation so that PjX = t{X)A for all X G Kp, where 
t{X) is a scalar that depends on X. After this adjust- 
ment, the stabilizer Kp is a codimension 1 subspace of 
the subspace of su(2)" spanned by Ai, A2, . . . , An- Thus 
there is some nonzero real vector (nii, m2, . . . , nin) such 
that 




tkAk : ^ mktk = 



In the following Proposition, we show that the only way 
for a nonproduct state to have this stabilizer is to be (LU 
equivalent to) a generalized n-qubit GHZ state. 

Proposition 6. Let p = be a pure nonproduct 

n-qubit density matrix with stabilizer subalgebra 



Then all the mj have 



.fe=i 



where ^ (toi, TO2, . . . , nin) G U 
the same absolute value and 10) is LU equivalent to a 
generalized n-qubit GHZ state. 
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Proof. Wc may suppose, after a suitable LU trans- 
formation, that Coo - o 7^ (the LU transformation 
exp(7r/2Cj) = Id (g) ■ • • ig) Id Q » Id ® • ■ • Id sends |/) to 
where Id is the 2x2 identity matrix and Ij denotes 
the multi-index obtained from / by complementing the 
jth index, and changes the stabilizer element ^ tkAk to 
'^{—iy^''tkAk, where Sjk denotes the Kronecker delta). 
This operation may change the sign of some of the m.j, 
but does not change their absolute values. 

If any mi = 0, then Ai € Kp. By Proposition [H it 
follows that the ith qubit is unentangled. But is not 
a product, so we can rule out the possibility that any me 
is zero. 

Suppose now that there exist two coordinates of 
(TOi,...,m„) with different absolute values. Then 
(toi, . . . , m„) and J ~ (ji, ■ • ■ , jn) are linearly indepen- 
dent, where J is any nonzero vector whose entries are all 
O's and I's. Let J be a multi- index not equal to 00 ■ • ■ 
with I's in positions ki, . . . , kp. We may choose a vec- 
tor (ti, . . . , t„) that is perpendicular to (mi, . . . , to„) but 
not perpendicular to J. It follows that X = J^^kAk lies 
in Kp, but X^fci ^fct 7^ Oi '^f' conclude that cj — 
by Proposition [21 This means that \tjj) is the completely 
unentangled state |00 • ■ • 0). But this is a contradiction, 
since \tjj) is not a product. This establishes that all the 
mj have the same absolute value. 

Thus \ip) has (possibly after an LU transformation) 
the stabilizer of a generalized n-qubit GHZ state, and 
therefore |?/;) must be a generalized 7i-qubit GHZ state 
by Proposition [5l □ 



V. THE 4-QUBIT CASE 

In this section we consider condition (ii) of Theorem [T] 
It follows from Lemma 1 and Proposition 1 of [ll| that 
after an LU transformation, if necessary, we may take Kp 
to be 

/ 4 4 ^ \ 

Kp=(Y^Ak,Y.B,,Y.cA. 

\k=l k=l k=l I 

From this assumption, we derive a canonical form for 
unique representatives of each LU equivalence class. 

Proposition 7. (Classification of 4-qubit states 
witii su(2) stabilizer) Let p = |?/') (V'l be a pure 4-qubit 
density matrix, where {ip) ~ '^iCi and let 

\k=l k=l fc=l / 

Then Kp = V if and only if 

IV') = a(|0011) + 11100)) 
+ 6(|1001)-F 10110)) 
+ c(|1010)-h 10101)) 



for .some complex coefficients a, b, c satisfying a > 0, 
abc =/= and a + b + c = 0. Furthermore, this is a unique 
representative of the LU equivalence class for \ip) of this 
form. 

Proof. Suppose that Kp = V. 

First we consider the consequences of the element 
J2^k in Kp. By Proposition [31 for each /, J for which 
pi^j = cjCj ^ 0, we have 

c(/,j)= Yl (-1^ = 0- 

From this it follows that if c/ , cj 7^ and wc have to flip 
no zeros and ui ones to transform / into J, then uq = ni . 
If m is the number of indices which are zeros in both / 
and J, then the number of zeros in / is tiq + m, and the 
number of zeros in J is 711 -I- m = tiq + m. It follows that 
the number of zeros is constant for all multi-indices / for 
which c/ 7^ 0. 

Next we consider J^^k in Kp. Since exp(7r/2C) = C, 
we have 

p = exp(^/2 ^ C,) • p = ^ PI.. J. \I) ( J| , 

so 

c/c} = c/cc}c (3) 

for all /, J. Applying ^ to the pair /, /, we obtain 
|c/| = |c/c|. Applying ([3]) to the pair 1,1'^, we see that 
c/c}e is real, so it must be that c/ — cic for all /. Since 
the number of zeros is constant for all multi-indices for 
nonzero state coefficients, it follows that this number 
must be two. Thus we have so far that 1-0) must be 
of the form 

= a(|0011)-h 11100)) 
+ 5(|1001) + 10110)) 
+ c(|1010) + 10101)) 

for some a, 6, c. If a = 0, then Ai — A2 would be in Kp, 
violating our assumption. Similarly, neither b nor c can 
be zero. 

To see that a+5+c must equal zero, note that (^ Ck)- 
p = means that (X^C'fc) ■ IV") = *^I0) for some t. A 
straightforward calculation shows that 

where denotes the multi-index obtained from / by 
complementing the fcth index. The coefficient of |1000) 
on the right hand side is (a + b+c) |1000), but zero on the 
left hand side, soa-|-6-|-c = 0. Finally, we may take a to 
be positive by applying a phase adjustment, if necessary. 
Conversely, suppose that ip has the form 

1^/-) = a(|0011) + |1100)) 
+ 5(|1001) + 10110)) 
+ c(|1010) + 10101)) 
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for some a, b, c, abc ^ and a + & + c = 0. It is easy 
to check that Kp contains V . By our analysis of sta- 
bihzer structure in [ll|, we have dimJi'p < 6. We know 
dim Kp = 6 if and only if is (LU equivalent to) a prod- 
uct of two singlet pairs ([1Q|), but we can rule this out 
(it is easy to check that I"!/)) is not a product). Lemma 1 
in [ll| rules out the possibility that dim Kp could be 4 or 
5, so dim iiTp must equal 3 and we must have Kp = V . 

To show that |^) is the unique representative of its 
LU class with the form stated in the Proposition, we use 
local unitary invariants. First, we consider invariants of 
the form tr((tryip)^) where ^ is a subsystem consisting 
of a subset of the 4 qubits. From these trace invariants, 
a straightforward derivation produces invariants h^h, h 
in the table below. 



Qubit pair in A 


Invariant 


(1,2) 


h{^)^\a\\h\ 


(1,3) 


= |a||c| 


(1,4) 


- \b\\c\ 



It follows that the norms of the state coefficients for 
1^) are LU invariant (for example, we have \a\ = 
^y Ii{ip)l2{'4')/ hitp))- Since we are taking a to be pos- 
itive, and since a + b + c = 0, a simple geometric or 
algebraic argument shows that b and c are determined 
up to conjugate. 

To determine the imaginary parts of b and c, we use a 
class of invariants given by 

p::rA^)= E cn---ci^c*j,---c*jr. 

where , . . . , is an m-tuple of 4-qubit multi- indices 
I'' = (ij, . . . ,i^), and 

jk — / -k -k -k -k \ 
•J — l*l,*cr(l),*r(l)'*</.(l)J 

for 1 < fc < 771 |T2| . We consider the invariant for m = 

3'^^ = (321)'^ = (213)' ^'''^ ^ ^ (231)- 
The imaginary part of this invariant acting on -0 is 

-2Aa'^bib2{bl + bl + abi) 

where 61 , 62 are the real and imaginary parts of b, re- 
spectively. From this we see that if 61 7^ 0, then 62 is 



determined. If 5i = 0, then 62 is not completely deter- 
mined; there remains the ambiguity of conjugation for b 
and c. To settle this, let a(r, s, t) denote the state 

r(|0011) + |1100)) 
-f- s(|1001) + 10110)) 
+ i(|1010) + 10101)) 

for state coefficients r,s,t. The (not local unitary) op- 
eration that permutes qubits 3 and 4 takes a{a, b, c) to 
a(a, c, b). If 6 is pure imaginary, then c is not, so our in- 
variant ^ . tells us that a(a, c, b) is not LU equivalent 
to a{a, c* ,b*). It follows that a{a, b, c) is not LU equiva- 
lent to a(a, b* , c*), for if the LU operation gi®g2®gz®gA 
takes a(a, 6, c) to a(a,6*,c*), then the LU operation 
51 ® 32 ® ff4 ® 53 takes a(a, c, b) to the LU inequivalent 
state a(a, c*, b*). 

Having exhausted all cases, we see that |?/') is the 
unique state in its LU equivalence class of the form stated 
in the Proposition. □ 



VI. CONCLUSION 

We have described the stabilizer subalgebra structure 
and have a complete description of states with stabilizers 
of maximum dimension, both for product and nonprod- 
uct states. We have shown evidence that warrants further 
study of stabilizer structure. There are at least three nat- 
ural directions in which to further pursue this analysis: 
seek complete descriptions, perhaps a classification, of all 
possible stabilizer subalgebras; seek complete description 
and perhaps canonical LU forms in the spirit of Proposi- 
tion [7] for states that have those stabilizers; and extend 
these pursuits to mixed states. 
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